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The one-dimensional kinetic Ising model introduced by Glauber in 1963 [T] is a 
prototypical model for the relaxation of a system towards thermal equilibrium. The 
model evolves by non-conservative single-spin-flip dynamics, with the requirement 
that the flipping rates fulfill the constraint of detailed balance with respect to the 
ferromagnetic energy (or Hamiltonian) 




( 1 . 1 ) 


in dimensionless units, where cr„ = ±1 are classical spins. The detailed balance 
condition implies that the dynamics is reversible, and hence that its stationary state 
is an equilibrium state. 

It has recently been realized [1 H il that the most general single-spin-flip 
dynamics for the ferromagnetic Ising chain with nearest-neighbour influence and spin 
reversal symmetry is a natural extension of Glauber dynamics in two directions 
associated respectively to non-linearity and irreversibility. The model introduced 


and 


in [HElIl] is defined by the rate functio n (|2.3[ ) (or equivalently ( |2.7[ )), which are 
extensions of the Glauber rate functions (2.11 and (2.3), where the parameters 5 
e take arbitrary values in the triangular region depicted in figure The 


stationary measure of this general model is that of the equilibrium one-dimensional 
Ising model, i.e., it is given by the usual Boltzmann-Gibbs distribution associated to 


the ferromagnetic Hamiltonian (1.1). 


The aim of the present work is to pursue the study undertaken in [3], where 
an exact analysis of the dynamical properties of the model along the line 5 = 0 of 
figures and was performed. The present study is concerned with the properties of 
the model at a generic point of the parameter space represented in figures and 
with particular emphasis on special lines and points. 

The main focus will be on the dynamical behaviour of two-spin correlation 
functions for a system started in either a random initial state or in a thermalized 
initial state. 


• In the first case all spin configurations are equally probable, each spin tT„ taking 
the values ±1 with probability 1/2 independently of the others, hence the two- 
time correlation 

C„(0,t) = (ao(0)a„(t)) (1.2) 

describes the transient regime of relaxation of the system to stationarity. We 
shall more particularly consider the autocorrelation 

C(0,t) = Co(0,t) = (ao(0)ao(t)), (1.3) 

which provides a measure of the overlap of the spin configuration of the system 
at time t with its random initial state. 

• In the second case the system remains stationary during its evolution, hence the 
correlation 

C'„,stat(i) = (o-o(0)cr„(t)) (1.4) 

gives a measure of the fluctuations of the system at stationarity. We shall more 
particularly consider the autocorrelation 

Ggtat (t) = C'o,stat(t) = (cro(O)cro(i))- 


(1.5) 
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We shall also provide some insight in the behaviour of the equal-time correlation 

Cnit) = {(TQ{f)an{t)) , ( 1 . 6 ) 

for the system initially prepared in a random initial state and relaxing towards 
stationarity. At long times this correlation converges to 

Cn = (croCT„) = (tanh/3)l"l, (1.7) 

denoting by /3 = 1/T the inverse temperature, which is the well-known expression of 
the correlation of the one-dimensional equilibrium Ising model [5]. 

The setup of the paper is as follows. Sections and provide some prerequisite 
knowledge on the model investigated in the present work. In section we give a 
reminder of the dynamical rules defining the model mill- In section ^ we give a 
reminder of some exact results [3] on the dynamics of this model on the solvability 
line ((5 = 0). We emphasize the existence of an oscillatory relaxation regime beyond 
a temperature-dependent threshold. Section is the bulk of the present work. It 
is devoted to novel aspects of the inhnite-temperature dynamics of the model, be it 
reversible or irreversible, beyond the solvability line. To this end we use a wide range of 
methods coming from various branches of statistical physics, each of these approaches 
shedding some light onto the problem from another angle. More specifically, we first 
use two general approaches, time series expansions (section |4.I[ ) and mapping of the 
dynamics onto a quantum spin chain (section |4.2[ ), before we analyze some special 
lines and points, including in particular the solvability line (i5 = 0) (section |4.3|), the 


reversibility line (e = 0) (section 4.4), the SEP point ((5 = —1, e = 0) (section 4.5), 


the dual SEP point ((5 = I, e = 0) (section 4.6) and the ASEP (microcanonical) 
line (i5 = —I) (section |4.7|). We also present some observations on the dynamical 


behaviour at a generic point (section 4.8), based on numerical simulations. Finally, 


an investigation of the spectra of the Markov matrix (section 4.9) provides a useful 
alternative tool to understand the qualitative features of the dynamics. In section 
we investigate the main novel features of the finite-temperature dynamics which were 
absent in the infinite-temperature situation. This includes a study of the relaxation 
rate along the reversibility line (section [5.1[ ) and the existence of an oscillating regime 
of relaxation beyond a threshold in the generic irreversible case (section [5.2| ). We then 
study the dynamics of the energy density of the model (section |5.3[ ) , and close with 
an investigation of two special points on the reversibility line (KDH and Metropolis) 
(section [A4| . We conclude by a brief discussion of our results in section]^ 


2. A reminder on the rate function of the model 


In this section we give an overview of the dynamical rules of the model introduced 
in 


2.1. The rate function of the model 

We first recall that the most general rate function for single-spin-flip dynamics 
((T„ —>■ — cr„) in continuous time with spatially homogeneous nearest-neighbour 
influence obeying both detailed balance and spin reversal symmetry reads [T] 

= -(^1 - -(1-H (5)cr„(cr„_i -h CT„+i) -f (5cr„_icr„+ij (Glauber). (2.1) 

X We refer the reader to those references for more details, in order to keep the reminder contained 
in this section succinct. 
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This rate function depends on two free parameters: a (which fixes the time scale) 
and 6. The parameter 7 is related to inverse temperature /3 by 

7 = tanh2/3, (2.2) 

as a consequence of detailed balance [1]. The so-called Glauber model usually refers to 
the case where the parameter 6 is set to zero. The success of Glauber model relies on its 
solvability, in the sense that its time-dependent behaviour can be determined exactly. 
All spin correlation functions of interest indeed obey linear evolution equations, that 
can be solved by analytical means. Whenever 5 is non zero, the solvability of the model 
is lost, in the sense that the hierarchy of evolution equations for spin correlations does 
not close. 


It has recently been shown mil] that the generalization of ( |2.1[ ) to a rate function 
only satisfying the condition of global balance reads 

= ^(1 -I- ecr„_iCT„ - ( 7(1 -I- 5) -I- e)tT„(j„+i -I- (Generic). (2.3) 

This rate function now depends on the additional parameter e. From now on, we set the 
time unit by fixing a = 1. Global balance is a weaker condition than detailed balance. 
It only requests that the (unique) stationary state of the system is given by the 
same Boltzmann-Gibbs measure as at equilibrium, corresponding to the ferromagnetic 


Hamiltonian (l.I) at fixed temperature T. A dynamics only obeying global balance is 


generically irreversible, and therefore generically leads to a nonequilibrium stationary 
state. Global balance contains detailed balance (leading to an equilibrium state) as a 
particular case. 

In the present context, the condition of detailed balance corresponds to fixing the 
parameter e to the value 


er--^(l- 


5 ) 


(where r stands for reversible), which yields ( 2.11 back, now rewritten as 
Wn = + Cr cr„((T„_i -|- cr„+i) -|- S tT„_icr„+i^ (Glauber). 

Whenever the condition e = Cr is not satisfied, the difference 


r] = e-e. 


(2.4) 


(2.5) 


( 2 . 6 ) 


quantifies both the irreversibility of the dynamics and its left-right asymmetry, as can 
be seen on the new form of (|2.3|) now rewritten as 


w, 


= ^ (1 + (Cr + r])crn-i(Jn + (cr “ ??)cr„cr„+i -|- 6 cr„_icr„+i^ (Generic) 


(2.7) 


which is invariant under the simultaneous change of ry into —77 and exchange of left 
and right. This property only holds for the Ising chain and does not extend to higher¬ 


dimensional situations [Ij. It turns out that the expression (2.3) (or equivalently (2.7)) 
not only satisfies the requirement of global balance once temperature is known, 
but represents the most general rate function for single-spin-flip dynamics with spin 


reversal symmetry and nearest-neighbour influence. Equation (2.2) now becomes a 
condition fixing the temperature [Ij. 








5 


Single-spin-flip dynamics of the Ising chain 


2.2. Spin-flip moves 


The spin-flip moves associated with each environment are listed in table The 
expressions of the rates corresponding to these moves in terms of the parameters 5 
and T] (arbitrary) and Cr (given by (2.4|) are deduced from (2.7). In the rightest 
column, the subscripts in the rates Wa-a' refer to the two neighbours a and cr' of a 

flipping (-I-) spin. The rates and w _(lines 1 and 4), which depend on S and 

Cl-, respectively correspond to the creation and the annihilation of a pair of domain 

walls, with energy cost AS = ±4. The rates w.^ _and w _(lines 2 and 3), which 

depend on 6 and g, correspond to the motion of the domain walls, respectively to the 
left and to the right, with no cost in energy (AS = 0). In other words, the parameter 
g = (w^ _— w _^)/2 is only involved in the asymmetric motion of the domain walls. 

An alternative description of the Ising chain can be given in terms of the domain 
walls [5] . Defining occupation numbers living on the bonds of the chain, according to 


= ^(1 - cr„tT„+i) = 0 or 1, 


( 2 . 8 ) 


then, if T„ = 1 (i.e., (T„cr„+i = —1), there is a particle (domain wall) on the 
corresponding bond; if t„ = 0 (i.e., cr„(T„+i = -1-1), there is a hole (i.e., no 
domain wall). On a finite system of N sites (and N bonds), with periodic boundary 
conditions, the total ferromagnetic energy £ is related to the number M of particles 
by £" = —N -\- 2M. The reactions among particles corresponding to the four spin-flip 
moves appear in the fourth column of table 


# 

spin-flip move 

spin-flip move 

reaction 

rate 

1 

-f-k-f ^ -f--k 

-^ --k- 

00 ^ 11 

W++ = ^(1 -|- 5) -|- Cr 

2 

-kH-H- 

-^---k-f 

01 ^ 10 

W+- = ^il-S)+g 

3 

--k-f ^-h 

H-^ -f-k- 

10 ^ 01 


4 

-H-^- 

H-h ^ -f-k-f 

11 ^ 00 

W _ = ^{1 -\- S) — Cr 


Table 1. List of spin-flip moves, with corresponding reactions among particles, 
and expressions of the corresponding rates in terms of th e parameters S and tj. 
The parameter Cr is related to 5 and temperature by ( |2.4[ |. 


The parameters 5 and e appearing in (2.3) are constrained by the condition that 
all rates must be positive. The parameter S has to obey \6\ < 1, while r/ = e — Cp has 
to obey \g\ < ry^ax, with 


f?max = ^(1 - <5). (2.9) 

This defines a triangular region of the S-e plane, depicted in figure 00 Two lines 
on this plot are remarkable. The first one is the reversibility line (e = Cr, i.e., g = 0), 
which corresponds to all possible reversible dynamics. The second one is the solvability 
line ((5 = 0), along which the dynamics is solvable, in the sense that spin correlation 
functions obey a closed system of evolution equations. These lines intersect at point G, 
representing the Glauber model {6 = g — 0). A detailed analysis of the dynamics along 
the solvability line, beyond the Glauber model, was given in |3] and will be recalled in 
section In figure the parameter g represents the vertical coordinate of a generic 
point relatively to the reversibility line e = Ci-. A 6-g plot gives a symmetric and 
temperature-independent representation of the triangle of all possible dynamics, as 
depicted in figure 
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Figure 1. Region of the 5—e plane where the dynamics obeys global balance 
and all rates are positive, with its two remarkable lines, the reversibility and 
solvability lines. These lines intersect at point G, representing the Glauber model 
(^S = rj = 0). The irreversibility parameter rj is the vertical coordinate of a generic 
point relatively to the reversibility line e = Cr. On the solvability line (5 = 0), 
correlation functions exhibit an oscillatory temporal behaviour in the out er re d 
segments, i.e., for e > 0 or e < —7, or, equivalently, for 770 < \7]\ < 1/2 (see 
Black squares show the SEP point (5 = —l,e = 0) and the two TASEP points 
(<5 = —1, e = ±1) (After [4].) 


For a generic point in the triangle shown in figures and the rates obey the 
unique constraint 


W++ _ „_4/3 _ 1 - 7 
w _ 1 + 7 


( 2 . 10 ) 


which is one of the two conditions for detailed balance. The other one, 


w+-=w-+, ( 2 . 11 ) 

only holds when the dynamics is symmetric or reversible (77 = 0). The boundaries of 
the triangle correspond to the vanishing of one of the rates given in table The left 

vertical side of the triangle {S = —1) corresponds to u>++ = w _= 0, so that only 

the motions of domain walls are allowed. Dynamics along this line will be hereafter 
referred to as microcanonical, as they conserve the ferromagnetic energy. On the two 

oblique lines one of the two conditions _= 0 (lower side) and w_^ = 0 (upper 

side) holds, i.e., p = =F 7 max- The point of intersection of these two oblique lines, to 
the right of the triangles, corresponds to _= w_= 0, i.e., <5 = 1 and 77 = 0, or 
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Figure 2. Same as figure^ represented in the S-rj plane. 


e = —Er = — 7 . The only allowed moves are the creation and the annihilation of pairs of 

domain walls, with rates tt>++ and w _respectively equal to 1=F7- More generally, the 

constant-energy moves, i.e., the motions of domain walls, are more favoured (resp. less 
favoured) than in the Glauber model for 5 < 0 (resp. i5 > 0) [7]. 

3. A reminder of the dynamics on the solvability line 

The solvability line (S = 0) is the natural extension to irreversible dynamics of 
the finite-temperature Glauber model (S = rj = 0). This line is parametrized by 
the irreversibility parameter 77 = e — Cr, with Cr = — 7 / 2 . Although the dynamics 
is asymmetric and irreversible, it is still solvable [3], in the sense that correlation 
functions obey a closed system of linear evolution equations, that can be solved 
analytically, just as for the Glauber model [D ISl IS- The outcomes are however 
surprisingly non trivial, as we recall below. 

3.1. Random initial state 

In the case where the system is started in a random initial state, the two-time 
correlation Cn{0,t) obeys the equation 

dC'„( 0 ,t) 


dt 


Cn{0,t) - (Cr - r/)Cn+l(0,t) - (Cr -f p)C„-i (0, t) , (3.1) 
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with C„(0,0) = Sn,o- Introducing the Fourier transform 

C(g,0,t) = ^C„(0,t)e-'"«, (3.2) 

n 

we readily obtain 

C(g,0,t) =e-^(«)*, (3.3) 

with 


r 2 (( 7 ) = 1 + 2(er cosg — iTysing). (3.4) 

Let us first consider the case where \r]\ < |er|. We have 

C„(0,t)=e-‘(^^j^) ' /„(2ty^?^), (3.5) 

where the /„ are the modified Bessel functions. In particular the autocorrelation reads 

C(0, t) = e-‘/o ( 2 t Ve 2 _ ,y 2 j , (3,6) 

The above results were derived in [3], where the ratio 77 /er is interpreted as a velocity 
and denoted by V. They demonstrate the existence of a threshold for the irreversibility 
parameter 77 , namely r\ = ±770 (i-e., V = ± 1 ), with 

m = \^A = l- (3.7) 

As long as I 77 I < 770 , the correlation C(0,t) decreases monotonically and falls off 
exponentially, with the decay rate 

Qfl = 1 — 2 a /£2 — 772 = 1 _ 1/^2 _ ( 3 ,g^ 

This rate is minimal at the Glauber point (77 = 0), where it equals 


Q!i = l- 7 , (3.9) 

and it increases towards 1 at threshold (see figure]^. 

Right at the threshold values 77 = ±770 (respectively corresponding to e = 0 and 
e = — 7 , the two endpoints of the green segment on figures and [^ , the dynamics 
become totally asymmetric, in the sense that the influence on the flipping spin comes 
from only one of the neighbours [Mill [311], with respectively 

Wn = ^(1- 7 cr„_icr„) and ■u;„ = ^(1 - 7 Cr„(T„+i). (3.10) 

Note that such a totally asymmetric dynamics can only occur if i5 = 0, since otherwise 
both neighbours always have an influence on the flipping spin. For 77 = — 770 , the 
correlation C'„(0,t) vanishes for n < 0, while for ti > 0 we have |3], 

C„(0,t) = ?7(^. (3.11) 

nl 

Likewise, for t] = tjq, the same correlation vanishes for n > 0, while for n < 0 we have 

c„(0,t) = ?;MT, (3.12) 

(-ny. 


The autocorrelation at the threshold, C(0,t) = e“*, is independent of temperature. 

Beyond threshold (770 < I 77 I < 1/2), i.e., in the outer red segments of figures 
and[^ the correlation (7(0, t) is given by the analytic continuation of (3.6), 

(3.13) 


(7(0,t) = e *Jo(2t^/rf 


- e; 
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where Jq is the usual Bessel function. This correlation now exhibits a damped 
oscillatory behaviour, in the form of asymptotically periodic oscillations, multiplying 
an exponential decay with unit rate. It vanishes for the first time at 


— 


J 


where j « 2.404825 is the first positive zero of Jq, while subsequent zeros are 
asymptotically separated by the half-period 


A = 


(3.15) 


2^/r]^ - e? 

The time scales ti and A both diverge according to the same inverse-square-root law 
as the threshold is approached (|? 7 | — rjo). 

The occurrence of an oscillatory behaviour beyond the threshold (Iryj > rjo) 
can be given the following simple explanation. Before threshold (|? 7 | < ? 7 o), both 
coefficients — (cr -I- ry) and — (cr — g) which appear in ( |3.1| ) are positive. The 
latter equation is therefore similar to a discrete diffusion equation, and leads to 
monotonically decaying correlations. Right at threshold (|ry| = go), one of the 
coefficients vanishes, giving rise to a totally asymmetric dynamics and to the totally 
directed correlations (3.11) and (3.12). Finally, beyond threshold ([jyl > go), the 
coefficients have opposite signs, giving rise to competing effects. The interpretation 
of (3.1) as a discrete diffusion equation is lost. The net outcome of these competing 
terms is the occurrence of damped oscillations. 


3.2. Thermalized initial state 


The case where the system is started in a thermalized initial state yields a richer 
behaviour, with two successive thresholds. 

The correlation C'„_stat(0, t), which accounts for the fluctuations of the system in 
the stati onary state, still obeys (3.1), albeit with the initial condition Cn^statCO, 0) = 
(see (1.7)), with v = tanh/?, so that 7 = 2v/{l -i-v"^), hence 

5stat(g,0) = --^—y. (3.16) 


I — 2v cos q 


The linearity of the differential equation (3.1) ensures that 

5stat(g,i) = dtat(g,0)e-^(«)‘. 


(3.17) 


In other words, 


C'ra.stat(0, t) is given by the spatial convolution 

/II \ 

C„.,tat(0,t) = e-‘^^;l"-’"l Im(2ty^ 



(3.18) 


We have in particular 


Cstat)!) = e 


,l™l 


Frl - g 

krl +?y 


ml2 




— g"^ 


(3.19) 


This correlation falls off exponentially whenever the irreversibility parameter g is less 
than the above threshold {\g\ < go). The corresponding decay rate however takes two 
different values in the following regimes. 
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• Regime /, corresponding to a weak violation of the equilibrium fluctuation- 
dissipation theorem [3], takes place for < rjc, with 


Vc = ^\/l - 72 


Here, C'stat(0 exhibits the same decay rate as C{ 0 ,t) (see (3.8|), i.e., 

Qfi = 1 — 


(3.20) 

(3.21) 


• Regime If corresponding to a strong violation of the fluctuation-dissipation theo¬ 
rem [3], takes place iov pc < \'ri\ < Vo- Here, the decay rate of C'stat(i) reads 

a2 = (3.22) 

7 

We have 02 < oii. 


Both rates match at 77 = pci as well as their derivatives with respect to 77, i.e., 

01 = 02 = 1 - 7 ^, ^ = ^ = {v = Vc)- ( 3 . 23 ) 

dr] dr] 7 

Figure]^ shows the 77 dependence of the decay rates ai and 02- Temperature is chosen 
in such a way that 77^ assumes its maximal value Pc = 1/4. This occurs for 7 = 
where we have ai(77c) = a2ipc) = 1/2, 770 = l/(2-\/2) and 02(770) = l/-\/2. 



Figure 3. Relaxation rates 01,2 against rj for (7(0, t) and C'stat(^)- In region I, ai 
holds for both correlation functions (in blue). In region II, cai holds for (7(0,t) 
(upper curve in blue), while 0:2 holds for Cstat(^) (lower curve in green). The 
black dotted vertical line is located at the critical value 77 c, which marks a 
bifurcation between two regimes of violation of the equilibrium fluctuation- 
dissipation theorem. The red dotted vertical line is located at the threshold 
value 770 beyond which both correlation functions oscillate. Temperature is chosen 
such that 7 = 1 /\/2. 


The behaviour of the correlation functions (7(0, t) (random initial state) and 
Cstat{t) (thermalized initial state) can be summarised as a phase diagram in the 7-77 
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plane, shown in figure]^ In Regime I (|? 7 | < pc): both correlations fall off monotonically 
to zero and share the same decay rate a\. In Regime II (rye < \r]\ < rjo), the correlation 
functions still fall off monotonically to zero, albeit with different decay rates, namely ai 
for C(0,t) and a 2 < oii for C'stat(t)- Finally, in the regions marked OSC (|? 7 | > po), 
both correlation functions exhibit oscillations. 



Figure 4. Dynamical phase diagram at finite temperature in the y-ry plane 
for the correlation C(0, t) and C'stat(t) on the solvability line {5 = 0). Black 
curves: g — ±?yc (see ( |3.20J I). Red lines: g = ±?yo (see ( |3.7[ )). The threshold 
values occurring in figureare recovered for 7 = 1/V2l 


In contrast, the equal-time correlation C'„(t) does not depend at all on the 
irreversibility parameter g [3]. This is a symmetry of the solvability line, which is 
(weakly) violated for generic values of the parameter 6 (see section 5.3). 

We refer the reader to [3] for more details on the results mentioned above and for 
further investigations on two-time observables. A parallel study has been performed 
both for the spherical model El and for the two-dimensional Ising model m- 


4. Infinite temperature 


We now address new aspects of the dynamics of the model, beyond the solvability 
line. In this section, which is the bulk of the present work, we focus our attention 
onto the case of infinite-temperature dynamics. In spite of its simplicity, this situation 
already encompasses most of the novel dynamical features of the generic model. The 
rate (2.3) becomes 


^I “f 1 “f ; 


(4.1) 


where the two free parameters 6 and e respectively encode the non-linearity and the 
irreversibility of the dynamics. The parameter e is indeed identical to g at infinite 
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temperature, thus figuredegenerates into figure]^ It has to obey |e| < emax, with 

(4.2) 


= Vn 




The dynamics defined by (4.1) is generically irreversible and gives birth to a 
nonequilibrium stationary state, except on the line e = 0, where it is reversible and 
yields an equilibrium state. This infinite-temperature stationary state is the random 
state where all spin configurations are equally probable. In other words, a random 
initial state is already thermalized. The correlation functions (7(0, t) and (7stat(i) are 
therefore identical. 

For the infinite-temperature Glauber model (J = 0, e = 0), the spins flip with 
constant rate 1/2, and thus remain independent of each other in the course of time. 
The two parameters of the model, S and e, deform this dynamics in a non-trivial 
fashion. The very simple statics of the model indeed does not preclude the occurrence 
of interesting dynamical features, which are investigated hereafter by means of a 
variety of techniques. We first use two general approaches, time series expansions 
(section |4.1[) 


and mapping of the dynamics onto a quantum spin chain (section 4.2), 
before we analyze some special lines and points, including in particular the solvability 
line ((5 = 0) (section the reversibility line (e = 0) (section [4(4| , the SEP point 
(<5 = —1, e = 0) (section [T^, the dual SEP point (S = 1, 


= 0) (section |4.6[ ) 

and the ASEP (microcanonical) line ((5 = —1) (section |4.7| ). We end up with some 
observations on the dynamical behaviour at a generic point (section 4.8), based on 
numerical simulations, and with an investigation of the spectra of the Markov matrix 
of the dynamics (section [T^. 


f.L Time series expansion 

In order to apprehend the role of the parameters S and e, a first approach consists 
in expanding the correlation (7(0, t) as a power series in time t. This technique is 
presented in full detail in |13| . In the present context, it will prove useful in identifying 
various symmetries of the dynamics. 

Consider the correlation (tTo(O)PA(i)), where 

PA{t) = n (4.3) 

n^A 

is the product of the spins of an arbitrary finite set A of sites. This correlation is 
non-vanishing only if the size |A| of the set A is an odd integer. It obeys the linear 
differential equation 

^(cro(O)PA(O) = -2 ^(cro(0)Wri(t)i74(i))- (4-4) 

Tie A 

The simplest of these equations, corresponding to A = {0}, reads explicitly 

^(cro(O)cro(t)) = - 2{aoi0)wo{t)ao{t)) 
at 

= - (CTo(O)cro(O) - <5(cro(0)tT_i(t)cro(t)cri(f)) 

+ e(CTo(0)cri(t)) - e(cro(0)cr_i(t)). (4.5) 

By considering larger and larger sets, such as A = { — I}, {!}, {—1,0,1}, and so on, 
and taking averages over the random initial state, one can systematically derive the 
coefficients 


Ofe = (-l)''^(o-o(0)cro(i))t=o 


(4.6) 
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of the time series expansion of the correlation of interest: 

C(0,t) = ^afc^. (4.7) 

fc >0 

A similar expansion can be derived for any other local observable. 

Symbolic routines run on the computer allow one to derive up to 20 coefficients 
of the above expansion. The first few of them read 


Oo = 1, 
ai = 1 , 


For 
(see 
in 6 


02 — 1 + (5^ — 2 e^, 

03 = 1 + 56 ^ — 6 e ^, 

04 = 1 + 18^2 3^4 _ ^^2 + 8S^)e^ + 6e^, 

05 = 1 + 58^^ + - (20 + 64(5^ + AS^)e‘^ + 30e^ (4.8) 

the infinite-temperature Glauber model (S = e = 0), we have C(0,t) = e“‘ 
(4.13)) so all the coefficients o^ equal unity. Generically the Ok are polynomials 
and e, whose degree grows linearly with the order k. 


4.2. Mapping onto a quantum spin chain 


Another technique which is very commonly used in investigations of kinetic Ising 
models [H HI in E] consists in mapping the dynamics of a chain of classical spins 
onto the statics of a quantum chain of spin operators in the spin- 1/2 representation, 


QX ^ oy ^ cfz ^ 


(4.9) 


where (t“, cr)(, tr^, are Pauli matrices acting at site nj§] 

There are in general several ways of mapping either the dynamics of a classical 
spin chain or a reaction-diffusion system onto a quantum Hamiltonian H. The review 
by Schiitz m presents a systematic way of doing so. Following this route, we obtain 
in the present case 


H = 




S(. 


'n^n-i-1 


+ o-^cr; 


n^n+l) + ^^i^n^n+1 



(4.10) 


On the infinite-temperature reversibility line (e = 0), H coincides with the 
Hamiltonian of the XXZ (anisotropic Heisenberg) quantum spin chain |181 , with <5 

being the anisotropy parameter. This spin chain is known to be integrable in the 
usual sense (existence of an infinity of conservation laws, applicability of Bethe Ansatz 


techniques). This property will be exploited in section 4.4 

Whenever the dynamics is irreversible (e 7 ^ 0), (4.10) is the Hamiltonian of the 
asymmetric XXZ chain, investigated in m- The latter Hamiltonian is still integrable, 
albeit non-Hermitian. It has complex spectrum in general (see section 4.9). For 
6 = —1, i.e., along the ASEP line, ferromagnetic interactions become isotropic. The 
corresponding Hamiltonian [20l [21] provides the basis for investigations of the ASEP 
by Bethe Ansatz techniques (see [12] for a review). 


§ This notation for the Pauli matrices should not be confused with the notation cjn for the classical 
Ising spin sitting at site n. 
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4-3. Solvability line (6 = 0) 

In this section we mention briefly how the results concerning the solvability line, 
recalled in section]^ simplify at infinite temperature. Equation (3.1) becomes 
dC'„( 0 ,t) 


dt 


= -C„(0,t) + e(C„+i(0,t) - C„_i(0,t)), 


(4.11) 


with the initial condition 0) = Sn,o- We thus obtain, in Fourier space, 

C{q, t) = with 

r 2 (( 7 ) = 1 — 2 iesin( 7 , (4.12) 

and finally Cn{0,t) = e~*Jn{2et). In particular the autocorrelation reads C{0,t) = 
e“‘Jo(2et). In the reversible case, i.e., for the infinite-temperature Glauber model 
{6 = e = 0 ), we are facing a dynamics of independent spins with rate 1 / 2 , hence 

C{0,t)=e-f (4.13) 

As soon as the dynamics is irreversible (e ^ 0), the correlation C{0,t) exhibits an 
exponential decay with unit rate, modulated by asymptotically periodic oscillations. 
The dynamics is indeed always in its oscillatory regime, because the threshold rjo 
vanishes at infinite temperature. The first zero of (7(0, t) occurs at time 

j 


t-\ — 


2|e| 


(4.14) 


where j « 2.404825 is the first positive zero of Jq, while subsequent zeros are 
asymptotically separated by the half-period 


A TT 


(4.15) 


/./. Reversibility line (e = 0) 

The dynamics is reversible on the line e = 0, where the second detailed balance 
condition (2.11) holds. This infinite-temperature reversible model has the following 


peculiarity. The coefficients Ok listed in (4.8) appear to involve only even powers of 


the non-linearity parameter 5. This suggests that the correlation (7(0, t) is symmetric 
under the change 5 —>■ —5. This is indeed an exact dynamical symmetry, which can be 
demonstrated as follows. Define new spin variables (?„, obtained from cr„ by flipping 
every second pair of spins, according tc|il] 

= (-l)i-‘(-/2)a„. (4.16) 

This construction was already used by Nemeth [7]. Consider now a given spin flip, 
expressed both in the original variables tT„ and in the new variables ct„. If the cr„ flip 
is of type 1 or 4 (see table [^, the flip is of type 2 or 3, and vice versa. In the 
infinite-temperature reversible case, flips of types 1 and 4 have the rate (l-|-5)/2, while 
flips of types 2 and 3 have the rate (1 — 5)/2. This explains the observed symmetry. 
The latter is broken as soon as the dynamics is irreversible (e ^ 0). This is testified 


by the presence of the term 4(5^e^ in the expression (4.8) of 05 . 


The spectrum of relaxation times of the model can be, at least in principle, 
extracted from the corresponding quantum Hamiltonian H. In particular the 
relaxation rate ai, characterizing the exponential fall-off of the correlation 


(7(0, t) 




(4.17) 


Int(fE) denotes the integer part of x, i.e., the largest integer less than or equal to x. 
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coincides with the spectral gap of H in the thermodynamic limit and in the appropriate 
sector. For the reversible model, where H is the integrable Hamiltonian of the XXZ 
spin chain, the gap in the relevant magnetic sector can be read off from |23j . In our 
units, it reads 

ai = Vl - < 52 . (4.18) 


We have thus found an exact expression for the relaxation rate of the reversible model 
at infinite temperature. 

For the infinite-temperature Glauber model (5 = 0), we have ai = 1, in agreement 
with the exponential decay C(0,t) = e“* of the correlation at all times (see (4.13)). 
More interestingly, the above relaxation rate vanishes with a square-root singularity 
at the endpoints of the reversibility line (5 —>■ ±1). The properties of the dynamics at 
these points are addressed in sections [4^ and |T6| 


4-5. SEP point (S = —1, e = 0) 


At this point, the dynamics looks particularly simple, as we have 

W++ = w _= 0, ui-i_= w _I- = 1. (4.19) 

This dynamics was already considered by several authors [24l [25] . In terms of the 
spin variables, it conserves the total energy E of the ferromagnetic model. It may 
therefore be referred to as a microcanonical dynamics. In terms of the particles 
representing domain walls, the dynamics conserves the total number M of particles. It 
coincides with the dynamics of the symmetric exclusion process (SEP) [T71I2H] . where 
the reversible diffusive moves 01 e-)- 10 occur with unit rate. The dynamics at this 
point is entirely independent of temperature, in agreement with its microcanonical 
character. 

The SEP point is one of the endpoints of the reversibility line, where the 
relaxation rate ai vanishes (see (4.18)), pointing toward a sub-exponential decay of 
the correlation C{0,t). The time series for the latter quantity involves coefficients 
which are pure numbers: 


oo = 1, ui = 1, 02 = 2, 03 = 6 , 04 = 22, Os = 90, 

06 = 396, 07 = 1848, 08 = 9108, og = 47400, ... (4.20) 

These numbers are not listed in the OEIS E3- It would be interesting to give them a 
combinatorial interpretation. 

The spin correlation (7(0, t) can be recast in the language of the SEP as follows. 
The spin ctq flips each time a particle crosses the origin. Let Q{t) be the net number of 
particles which cross the origin during a lapse of time of duration t (counted positively 
if moving to the right, and negatively if moving to the left) in the stationary state of 
the SEP characterised by a particle density p (a random initial spin state corresponds 
to p = 1/2). We have 

C(0,t) = ((-1)'5(*)) = (e*"«(‘)). (4.21) 

After a long time t, Qft) will be typically large, and hence approximately given by 
Q{t) « pi?(t), where R{t) is the random position at time t of the particle which was 
the first to the right of the origin at time t = 0, say. The distribution of R{t) has 
been studied by several authors [2S]. The bulk of this distribution is known to be 
asymptotically a centered Gaussian, with a variance growing as 

/l. (4.22) 

P V TT 
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Taking this result literally, forgetting about the discrete nature of particles, we obtain 
the rough estimate 

(7(0,1) - exp ~ exp (-p{l - p)VnHj . (4.23) 

A stretched exponential decay of the correlation (7(0,1), of the form 


(7(0,1) ~ exp (^—A{p)Vi 


(4.24) 

has indeed been predicted in [^, where the amplitude A{p) is expressed in terms of 
the solution of a variational problem. A quantitative prediction for this amplitude, 

(4p(l - 




k>l 


fc3/2 


(4.25) 


can be read off from the work by Derrida and Gerschenfeld |53], involving the Bethe 
Ansatz and the use of results by Tracy and Widom. 

Coming back to the Ising chain, a random initial state corresponds to a particle 
density p = 1/2, where the amplitude A{p) takes its maximal value 


A(l/2) 


C( 3 / 2 ) 


1.473874, 


(4.26) 


where C is the Riemann zeta function. Interestingly enough, the particle density 
p = 1/2 appears as a singular point, around which the amplitude A[p) has a triangular 
shape: 

A(p)« A(l/2)-4|p-l/2|. (4.27) 


4-6. Dual SEP point (S = 1, e = 0) 


This point, where = w _= 1 and _= w _= 0, is the other endpoint of 

the infinite-temperature reversibility line. In terms of the spin variables (t„, all moves 
which change the ferromagnetic energy £, i.e., moves of types 1 and 4 (see table[^, are 
equally allowed and take place with unit rate. In terms of the particles representing 
domain walls, the dynamics consists in the pair creation and annihilation reactions 
00 (->■ 11 with unit rate. This dynamics has a simple alternative description. In terms 
of the new (dual) particles 


— 2 ^^ ^nf^n-t-l) — 


Tn 
1 - 


(n even) 
(n odd) 


(4.28) 


(see (4.16)), the dynamics is again that of a SEP. 

The symmetry related to the transformations (4.16) and (4.28) of spins and 
particles exchanges the SEP and dual SEP points. In particular the correlation (7(0, t) 
is identical at both points. We recall that this symmetry only holds in the infinite- 
temperature case. 
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4 . 7 . ASEP (microcanonical) line (5 = —\) 

Along this whole line, we have w++ = w _= 0, therefore the dynamics is 

microcanonical. In terms of the particles representing domain walls, the dynamics 
consists in the moves 01 —>■ 10 and 10 —>■ 01 with respective rates 1 + e and 1 — e. 
Particles go preferentially to the left for e > 0 and to the right for e < 0. This 
dynamics is that of the asymmetric exclusion process (ASEP). It is invariant under 
the simultaneous change of e into —e and exchange of left and right. The reversible 
dynamics of the SEP point is recovered when e = 0. The endpoints e = ±1 correspond 
to the totally asymmetric exclusion process (TASEP) [T71 [2^ . 

The time series of the correlation C{Q, t) at both TASEP points again involves 
coefficients which are positive integers. 


oo = 1, ai = 1, 02 =0, 03 = 0, 04 = 8 , 05 = 40, 

Og = 136, 07 = 392, og = 1032, 09 = 3912, ... 


(4.29) 


For arbitrary values of e, the correlation C{0,t) is still related by the iden¬ 
tity (4.21) to the net number Q{f) of particles which cross the origin during a lapse 
of time of duration t in the stationary state of the ASEP with a particle density p. 
A random initial spin state again corresponds to p = 1/2. The distribution of Q{t) 
is however no longer Gaussian, as was the case for the SEP. It is indeed governed 
by the non-linear Kardar-Parisi-Zhang (KPZ) theory [30] ■ The mean value of Q{t) 
is given by the mean stationary current, i.e., {Q{t)) = Jt, with J = —2p(l — p)e in 
our units. The magnitude of the fluctuations of Q{t) around its mean value depends 
on the way it is measured. For a given particle and a given configuration of the 
system at time < = 0, typical fluctuations scale as (e^t)^/^. Their statistics are by now 
well characterised m- they involve one of the universal laws originally discovered 
by Tracy and Widom in random matrix theory. The situation is however made more 
intricate by the fact that intrinsic fluctuations are masked by statistical ones, as soon 
as they are averaged either over different particles and/or over different configurations 
of the system at time t = 0 [33|. To the best of our knowledge, no analogue of the 


results (4.24), (4.25) has been derived for the ASEP so far. 


We have investigated the correlation (7(0, t) by means of numerical simulations. 
It is observed to exhibit a stretched exponential law of the form 


C(0,t) ^ exp , 


(4.30) 


i.e., the same functional form as the exact result (4.24) which holds at the symmetric 
SEP point (e = 0). The amplitude A depends continuously on the irreversibility 
parameter e. The decay law (4.30) is modulated by oscillations, whose half-period A 
also depends on e (see figure [6 below). The exponent 1/3 characterising the anomalous 
scaling of intrinsic fluctuations in one-dimensional KPZ theory does not enter (4.30). 


This feature, which may seem surprising at first sight, is certainly related to the above 
discussion on the nature of the fluctuations. It would be very desirable to obtain a 


derivation of our conjectured result (4.30) from the vast body of knowledge on KPZ 
theory. 

Let us first present our numerical data for the TASEP point (e = 1). Figure]^ 
shows absolute logarithmic plots of (7(0, t) against time t (left) and against \/i (right). 
The left panel shows the periodic pattern of oscillations, with a first zero at « 1.26 
and a half-period A « 1.84. The slope of the straight line drawn on the right panel 
yields A « 2.06 for e = 1. 
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Figure 5. Absolute logarithmic plots of the correlation (7(0, t) at the 
TASEP point. Left: plot against time t. Bar with left-right arrows: half¬ 
period A ~ 1.84. Right: plot against %/t. The straight line tangent to the 
data has an absolute slope A « 2.06. 


Repeating the same analysis for several positive values of e along the ASEP line, 
we obtain the data shown in figure The left panel shows the amplitude A of 
the stretched exponential law ( 4.30[ ) against e. This quantity is observed to depart 
continuously from the analytically known SEP value (see (4.26)), shown as a blue 
square. The right panel shows the products eti and eA against e, where ti is the 
location of the first zero of (7(0, t), while the semi-period A is the asymptotic lapse 
of time between two consecutive zeros. Both products hardly depend on e. So, the 
scaling laws 

(4.31) 

e 

which hold exactly along the solvability line (see ( |4.14 1, (4.15)), also hold approxima¬ 
tely along the ASEP line. 



Figure 6. Left: amplitude A of the stretched exponential law (4.301 along 
the ASEP line, against e. Blue square: analytically known SEP value at 
e = 0 (see ( |4.26[ )). Right: products eti and eA along the ASEP line, 
against e. 
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We now turn to a brief description of the infinite-temperature dynamics for generic 
values of the parameters 5 and e inside the triangular region shown in figures and 
The correlation (7(0, t) falls off exponentially, as 

(7(0,t)-e-“i‘, (4.32) 


where the decay rate a,\ has a rather weak dependence on the irreversibility 
parameter e. This exponential decay is modulated by oscillations which are 
asymptotically periodic in time. Here again, both the first zero ti and the half¬ 
period A roughly follow the 1 /e law ( |4.31[ ). 

These features are illustrated in figure for S = —1/2, hence Cmax = 3/4. Data 
are only given whenever the plotted quantities can be measured with a reasonable 
enough accuracy, i.e., for not too small values of the irreversibility parameter e. Even 
so, statistical errors are more important than along the ASEP line (compare figures]^ 
an d^. The left panel shows the decay rate ai against e. This quantity s eems t o depart 
continuously from the analytically known value ai = at e = 0 (see (4.18)), shown 
as a blue square. The right panel shows the products eti and eA against e. These 
products again exhibit a very weak dependence on e, so the scaling laws (4.31) hold 
approximately for generic values of the parameters b and e. 



Figure 7. Left: decay rate ai again st e for & = —1/2. Blue square: known 
value ai = y/2>l2 at e = 0 (see (4.181). Right: products eti and eA against e 
for 5 = —1/2. 


4-9. Spectrum of the Markov matrix 

A useful alternative way of investigating the problem consists in looking at the 
spectrum of the Markov matrix which represents the generator of the stochastic 
dynamics on a finite chain. This approach has already proved to be a useful tool 
in several circumstances [331 IMl ISS]. Here again, the shape of the spectra of Markov 
matrices will provide a clear picture of the main characteristics of the dynamics, 
including reversibility or integrability. 

A finite system of N sites with periodic boundary conditions has 2^ configurations 
{cti ... ctat}, that we assume to be ordered lexicographically from {-!- -I- -I-• • •} to 
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{-- •}. Let P(t) denote the vector of the probabilities of these configurations at 

time t. This vector obeys the differential equation 




(4.33) 


where the 2^ x 2^ matrix M is the Markov matrix of the problem. Its non-diagonal 
entries are the transition rates between adjacent configurations (i.e., configurations 
which differ by a single spin flip), while its diagonal entries are determined by the rule 
that column sums vanish, ensuring the conservation of probability. 

If the dynamics is reversible (e = 0, i.e., _= w _|-), the Markov matrix M 

is symmetric and its eigenvalues are real. Otherwise M has complex spectrum in 
general. Let us henceforth denote the eigenvalues of M as — A^, with Re > 0. As the 
dynamics obeys spin reversal symmetry, the eigenvectors of M have a definite parity, 
i.e., they are either even or odd under spin reversal. Even and odd subspaces have 
equal dimensions 2^“^. The eigenvalue Ai = 0, corresponding to the stationary state, 
pertains to the even sector. The stationary state is unique, so Ai is non-degenerate, 
except in the extremal cases (<5 = ±1) which obey conservation laws. Finally, the 
spectra of the Markov matrix M and of the corresponding quantum Hamiltonian H 
(see section 4.2) are related as follows: the spectrum of H consists of two copies of 


the spectrum of M in the even sector, while the eigenvalues of M in the odd sector 
do not appear in H. 

The smallest generic system consists of = 3 sites. The corresponding 8x8 
Markov matrix M and its eigenvalues can be written down explicitly. We have 

\ 


M = 


-3w++ 

W++ 

W++ 

0 

W++ 

0 

0 

0 

W++ 

—s 

0 

W- + 

0 

w+- 

0 

0 

W++ 

0 

—s 

w+- 

0 

0 

W- + 

0 

0 

w;+_ 

W- + 

—s 

0 

0 

0 

u;++ 

W++ 

0 

0 

0 

—s 

W- + 

W+- 

0 

0 

W- + 

0 

0 

w+- 

—s 

0 

W++ 

0 

0 

W+- 

0 

W-+ 

0 

—s 

W++ 

0 

0 

0 

W++ 

0 

W++ 

W++ 

-iw++ 


(4.34) 


where the rates rc++, _, w _^ are given in table in the infinite-temperature case 

(er = 0), and with the shorthand s = rc++ -f _-I- w _The eigenvalues of M are 

as follows. 

Even sector: Ai = 0, A 2 = 2(1-|-15), A 3_4 = 2 — J ± i e-\/3. 

Odd sector: As^g = 2 ± \/l -|- 3(5^, Ay^g = 1 ± ieVZ. (4.35) 

These expressions are not invariant under the change 5 —> —5. This symmetry of the 
infinite chain indeed only holds on systems whose size iV is a multiple of 4. 

Let us now describe the spectrum of the Markov matrix of infinite-temperature 
dynamics on larger systems. The following discussion (see the contrast between 
figures andclearly demonstrates that an investigation of these spectra provides 
a very sensitive tool to detect characteristic features of the underlying dynamics, such 
as integrability or irreversibility. 

• Infinite-temperature Glauber model (6 = e = 0). As already mentioned, this case 
corresponds to a dynamics of independent spins. The result (14.131) generalizes 


as follows. The correlation function built on any fcuple of distinct spins decays 
exponentially as 

(cr,i (t)... CTi, {t)ai, (0)... CTi, (0)) = (4.36) 







Im A 
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The spectrum of the Markov matrix is therefore highly degenerate |35] : it consists 
of the integers k = 0,..., N, with the combinatorial multiplicities 


fi{N, k) 



N\ 

k\[N -k)V 


(4.37) 


The even (resp. odd) sector corresponds to even (resp. odd) values of k. 


• Solvability line (5 = Q). Along this line, parametrized by the irreversibility 
parameter e, the eigenvalues of the Markov matrix are still given by simple 
formulas. They are indeed observed to be linear combinations, with integer 
coefficients, of the following complex frequencies 

71 TT 

0„ = l-il„, /„ = 2esin—. (4.38) 


The latter frequencies are obtained by inserting into the dispersion law (4.12) 
discrete quantized momenta of the form q = mr/N, corresponding to both 
periodic and anti-periodic boundary conditions. So, the real part of any eigenvalue 
is still an integer k = 0,..., N, while its imaginary part is an integer linear 
combination of the Z„. It is therefore strictly proportional to e. Figure shows 
the upper parts (Im A > 0) of the spectra of the Markov matrices for N = 6 and 
N = 7, with 5 = 0 and e = emax = 1/2. The imaginary parts of the eigenvalues, 
shown by horizontal lines, read in increasing order (bottom to top): For N = 6 
(left): h, I 2 , I 3 , h + h, For N = 7 (right): I 3 — I 2 , h h — h, h — h, h, 
h — h, h + I 3 — h, ^ 2 , I 3 , h + h, h + I 3 — h, h + h, h + h, h + h + h- 


“8 - 1 - i - ?- 1 - ?- 

Re A 



Figure 8. Upper parts (Im A > 0) of the spectra of the infinite-temperature 
Markov matrices M for A = 6 (left) and A = 7 (right) with 5 = 0 and 
e = Cmax = 1/2. Blue symbols: even sector. Red symbols: odd sector. 
Horizontal lines: imaginary parts listed in the text. 


• Reversibility line (e = 0/. Along this line, parametrized by the non-linearity 
parameter 6, the spectrum of the Markov matrix is real. Figure shows this 
spectrum for A = 8 against S. The plot is invariant under the change 5 —>■ — 5, 
as expected as the system size is a multiple of 4. The highly degenerate integer 
spectrum of the Glauber point is manifest in the middle of the plot (5 = 0). Level 
crossings take place in both sectors. The lowest eigenvalue A = 0, corresponding 
to the stationary state, is non-degenerate, except at the endpoints (5 = ±1), 
where the conservation laws induce degeneracies. 
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Figure 9. Spectrum of the Markov matrix M of the infinite-temperature 
reversible dynamics against the non-linearity parameter 5 for N = 8. Blue 
lines; even sector. Red lines: odd sector. 


• Generic case. For arbitrary parameter values ((5 7 ^ 0 and e 7 ^ 0), the spectrum 
of the Markov matrix appears as a rather structureless cloud in the complex A- 
plane. Figure shows this spectrum for iV = 12 in two cases. At the TASEP 
point (5 = —1, e = 1) (left), eigenvalues tend to accumulate near ReA = 0. 
This observation goes hand in hand with the fact that the eigenvalue A = 0 is 
degenerate, because the dynamics conserves the total energy. At a generic point 
{6 = —1/2, e = 2/3) (right), there are fewer eigenvalues near ReA = 0. The 
outermost part of the spectrum shows alternating blue and red stripes, which are 
remnants of the ordered structure of the spectrum in the solvable case. 



8 10 12 14 16 18 0 2 4 6 8 

Re A Re A 


Figure 10. Spectrum in the complex A-plane of the infinite-temperature 
Markov matrix M for N = 12. Left: TASEP point. Right: a generic point 
{S = —1/2, e = 2/3). Blue symbols: even sector. Red symbols: odd sector. 
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In this section we focus our attention onto the novel features of the finite-temperature 
dynamics which were absent in the infinite-temperature situation, investigated in 
section |4l 


5.1. Reversibility line (r] = 0) 

The finite-temperature reversible dynamics exhibits several novel features with respect 


to its infinite-temperature counterpart (see section 4.4). 

The correlation functions C(0,t) (random initial state) and C'stat(t) (thermalized 
initial state) are different from each other. The analysis of the solvable case {S = 0) 
recalled in section however strongly suggests that, for a reversible dynamics, both 
correlations fall off exponentially with a common decay rate ai, which depends both 
on temperature and on the non-linearity parameter S. 


The first few coefficients of the time series expansion (4.7) of the correlation 
function (7(0, t) read 

Oq = 1 , 
oi = 1, 

02 = 1 T 2e^ -|- (5^, 

03 = 1 -l- 6 e^ -|- 8e^i5 -I- 55^, 


04 = 1 + 126 ? + 6 ef -h 48e?(5 -h 2(9 -h 16e?)(5^ -f 3(5'‘ 


(5.1) 


The invariance of the infinite-temperature reversible dynamics under the change 
6 —> —6 is broken at any finite temperature. This is testified by the presence of a 


term proportional to e^S in the coefficient 03 in (5.1). Therefore, at variance with 


the infinite-temperature situation (see (4.18)), the decay rate ai is not expected to 


be even in S. In particular the dynamics at the two endpoints (<5 = ±1) are now 
different. For S = —1, we recover the SEP point and its microcanonical dynamics, 
investigated in section [431 irrespective of temperature. For <5=1, the allowed moves 
are the creation and annihilation of pairs of domain walls, with respective rates 1 T 7 - 
This temperature-dependent dynamics conserves the total number of dual particles 


described by the occupation numbers t„ (see (4.28)). In the zero-temperature limit, 
this non-generic dynamics becomes an interesting example of a kinetically constrained 
model, which exhibits all the generic features of metastability. This model has been 
investigated in detail in |36j . together with various other one-dimensional examples. 

The corresponding quantum Hamiltonian reads 


H = 


2 ^ 


. v^S — 1 
1 + ^ 0 - 7 ct: 


1 + v^ 


n^n+1 


-k 


5 — 


\+v 


+5ala. 


2 ^n^n+1 


n^n+1 


2u(l -k5) 
1 + v'^ 


(5.2) 


where we recall that u is a shorthand notation for tanh /3. This Hamiltonian was 
already derived by Nemeth [7] (up to a global factor 2 and with a different convention 
for the sign of erf). It describes the XYZ (fully anisotropic Heisenberg) spin chain 
in a uniform external field. This model is known not to be integrable in general. As 
a consequence, we have no analytical prediction for the relaxation rate ai, except 
in the solvable case {6 = 0), where a\ = 1 — 7 (see (3.9)), and at the endpoints 


(<5 = ±1), where ferromagnetic interactions become isotropic, so that H becomes 
invariant under spin rotations around the 2 : axis. This symmetry goes hand in hand 
with the conservation laws. As already underlined in [7], it implies the vanishing of ai 
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at these points. Furthermore, the relaxation rate ai can also be determined exactly 
at a special point known as the KDH point (see section 5.4). Finally, yet other 
special cases of the Hamiltonian (5.2) are integrable as they correspond to physical 
realizations of Hecke algebras m- The latter cases however do not bring more results 
on the present problem. 

In order to investigate the dependence of the relaxation rate ai on the non¬ 
linearity parameter 6 , we have measured the correlation (7(0, t) by means of numerical 
simulations. Figure [IT] shows our results for a case of moderate temperature (v = 1/2, 
i.e., 7 = 4/5). The plotted values of ax are reasonably accurate. The estimated error 
bar is comparable to the symbol size on the figure. Furthermore, we have checked 
that the stationary correlation Cstatif) yields compatible values of ai within the error 
bar. The extrapolation procedure looses its accuracy near the endpoints (S —>■ ±1), 
where ai is known to vanish. This is to be expected, as the correlation exhibits a 
crossover to the stretched exponential law (4.24) as (5 —> — 1, and to a similar kind 
of subexponential relaxation law as (5 —>■ 1. The plotted values for ai seem to vanish 
according to the same square-root law as the infinite-temperature result (4.18). They 
are also compatible with the presence of a slight cusp at the solvable point ((5 = 0), 
shown as a blue symbol, where ai = l — 7 = 1/5. The full curves show the outcomes 
of two separate fits of the data points for 5 < 0 and (5 > 0, where ai is assumed to be 
the product of the result (4.18) by a second-degree polynomial. 



Figure 11 . Relaxation rate oi of the correlations (7(0, t) and (7stat(t) on 
the finite-temperature reversibility line fq = 0 ) against the non-linearity 
parameter 5 for v = 1/2, i.e., 7 = 4/5. Blue symbol: exact result 
ai = 1 — 7 = 1/5 in the solvable case (S = 0). Green symbol: exact 
result for the KDH point (see ( |5.19[ )). Red symbols: numerical data. Full 
curves: fits (see text). 
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In the generic situation of an irreversible dynamics, the most natural question is 
whether the non-trivial phase diagram in the 5-rj plane depicted in figure subsists 
at finite temperature. 

We have evidenced, by means of numerical simulations, the existence of the 
threshold ? 7 o between a regime where the correlation (7(0, t) falls off monotonically 
and a regime where it is oscillatory. We have not addressed the more delicate question 
(from a mere numerical standpoint) of the existence of the other threshold ijc, related 
to the non-analyticity of the decay rate of the thermalized correlation Csta,t{t). 

The threshold po between the oscillatory and monotonic regimes has been 
determined as the value of p where ti, the first zero of the correlation (7(0, t), diverges. 
Figure 12 shows l/t\ against p^ for 7 = 1/2 and 5 = —0.7. This way of analyzing 
data is inspired from the analytical result (3.14) in the solvable case (<5 = 0). The 
data points (red symbols) exhibit an almost perfect linear law. A quadratic fit (full 
curve) yields the threshold po = 0.162. 



Figure 12 . Plot of 1/tf against for 7 = 1/2 and <5 = —0.7. Red symbols: 
data points obtained by numerical simulation. Full curve: quadratic fit 
yielding the threshold 70 = 0.162. 


Repeating the same analysis for several values of <5, we obtain the phase diagram 
in the S-p plane shown in figure 13 for 7 = 1 / 2 . Symbols show the dependence 


on S of our prediction for the threshold po (and its symmetric counterpart —po). The 
oscillatory regions {\p\ > po) are marked OSC. The threshold po is observed to be 
nearly constant for d > 0 and equal to its value r 7 o= 7/2 = l/4in the solvable case 
(d = 0) (see (3.7|), shown as green squares. The ht (full curve) suggests that the 
threshold vanishes with a square-root singularity. 


Po ~ Vl S, 

as the SEP point is approached (<5 —>■ —1). 


(5.3) 
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Figure 13. Phase diagram in the 5-g plane for 7 = 1/2. Red symbols: data 
for the threshold ±70• Full black curves: fits suggesting the square-root 
law (5.31. Green squares: exact values ±r;o = ±1/4 for <5 = 0 (see figures]^ 
and|^. 


To close, let us underline that the existence of the threshold tjq in the d-rj plane 
is not visible on the spectra of Markov matrices, which always look like the infinite- 
temperature spectra shown in figure [TOl 


5.3. Equal-time spin-correlation function 


Another novel feature of finite-temperature dynamics is the non-triviality of the equal¬ 
time spin-correlation function Cn{t) (see (1.6)), starting from a random initial state. 

Here we focus our attention onto the correlation between neighboring spins 
C'i(t) = (o'o(t)o'i(t)), i.e., (minus) the ferromagnetic energy density per site (or bond). 
At infinite temperature, for an arbitrary spatially homogeneous initial state, this 
quantity obeys the closed evolution equation 


dCijt) 

dt 


-2{l + S)Ci{t). 


(5.4) 


For a random initial condition ((71(0) =0), we thus recover that Ci{t) remains equal 
to zero all throughout its evolution. 

At finite temperature, in contrast, Ci(t) increases monotonically from zero 
towards its equilibrium value Ci = tanh,d (see O)- Its time series expansion reads 

(7i(t) = 5]6fc(-l)'=-i^, (5.5) 

A :>1 


with 


bi =7(1±(5), 

62 = 27(l±5)^ 
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h =7(1+ <5)^(4+ 7^), 

64 = 27(1 + 5 )^( 4(1 + 7^)(1 + ( 5 ) - 7^), 


(5.6) 


It is striking to observe that there is no dependence in the irreversibility 
parameter 77 before the 7th order. At this order, 67 contains a term of the form 
87^(1 + Similarly, bs contains a term of the form 247^(1 + <5)^(8 + 5)5'^rf, 

while 69 contains one term in rY and one term in rf. As it turns out, rY is always 
accompanied by 5'^, as the correlation Cn{t) does not depend on the irreversibility 
parameter rj at all in the solvable case {5 = 0) (see section 3.2). A similar very weak 
effect of the irreversibility on equal-time spin correlations has been observed both 
numerically and by short-time expansions for the two-dimensional Ising model in its 
high-temperature phase m- 

At high temperature, the full temporal behaviour of C'i(t) can be derived as 


follows. To first order in 7 , the coefficients of the time series expansion (5.5) assume 
the simple form ~ 2^“^(1 -|- SY'y, irrespective of rj. This has been checked up to 
the 9th order. Hence 


Ciit) « ^ . 


(5.7) 


The asymptotic value Ci « 7/2 agrees with the exact result, i.e., u, as 7 « 2u « 2/3 
to leading order at high temperature. The above result suggests that the convergence 
rate of the energy, such that 


tanh /3 — Cl (t) ~ e 


-cke* 


(5.8) 


takes the simple value 


aE = 2{l + S) (5.9) 

to leading order in the high-temperature regime. This result coincides with the 


prefactor in the right-hand side of (5.4). It vanishes when d = — 1, as should be, 


since this corresponds to the SEP point whose microcanonical dynamics conserves the 
total energy. Finally, the rate a_E is also mentioned in [23] as being one of the gaps 
governing the low-temperature thermodynamics of the quantum XXZ model. 

At finite temperature, we have no analytical prediction for aE in general, which 
is expected to have a (weak) dependence on the irreversibility parameter rj. In the 
solvable case {6 = 0 ), we have as = 2 ai = 2(1 — 7 ), irrespective of p l^. 


5.4- A special point on the reversibility line 

The dynamics simplifies for the special point considered by Kimball 
and Haake (33, where 


5 = 


7 


2 - 7 ’ 


yielding 


1 / ^ 

~ 1 (^n(<^n—1 ^n+l) 

2 \ 2--f 


and by Deker 

(5.10) 

(5.11) 


In the S-e plane, this point is located at the intersection of the reversibility line 
with the second bisectrix, and therefore obeys the condition 6 = —Ci-. The KDH 
point is alternatively characterised by the fact that three of the four rates coincide: 


W++ = W-i_= w _as can be seen on (5.11). 
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For an arbitrary reversible dynamics, and an arbitrary magnetized initial state, 
the evolution equation for the magnetization profile (cr„) reads 


= -(cr„) - er((cr„-l) + (ct„+i)) - S {-Kn), 


where 


'^n — ^n—l^n^n+1- 

The evolution equation of the latter quantity reads 

= -3(7r„) - 2er(((T„_i) + {(7n+l)) “ <5 (ct„) + {(j)n), 


(5.12) 

(5.13) 

(5.14) 


where 

4^n — €r(f^n—2f^n<^n+l T —l^n<^n+2) l^n+lf^n+2 “t" ^n—2^n—l^n-\-l)- (5.15) 

Consider a system of N sites with periodic boundary conditions. If Cr = —S, the 
spatial sum of the (pn vanishes identically. Defining 

n(t) = ^ 

n n 

we obtain from (5.12) and (|5.14[) two coupled linear equations for these quantities: 


M{t)\ _ (25-1 -5\(M{t)\ 
dt (vH)!) )~\i5 " ■ I 


(5.17) 


-3An(t);- 

The corresponding decay rates are the opposites of the eigenvalues of the above matrix, 

i.e.EHi™, 


a± 


= 2 - (5 ± a/I + 2d - 2(52 = 


4 - 37 ± ^4 - 372 

^ 2 A ■ 


(5.18) 


For a ferromagnetic model, it is physically reasonable to assume that the spin 
autocorrelation (7(0,1) falls off at the same rate as the total magnetization M{f). 
This line of thought yields the identity 

Qfi = a-. (5.19) 

For 7 = 4/5, the KDH point is at (5 = 2/3 and (5.19) reads Oi = (4 — •\/l3)/3 = 
0.131482. This prediction, shown in figure E] as a green symbol, agrees with our 
numerical results. 

It has been noticed in SOI that a similar construction works for the reversible 
dynamics with 

(5.20) 


5 = - 

2 + 7 

Let us point out that this choice corresponds to the Metropolis dynamics: 
Wn = min(l,e“^^^) 


2 + 7 
2(1 + 7) 


^1 2 -)- y + f^n+l) + ■ 


(5.21) 


Here again, three of the four rates coincide: _= w _^ = w _The Metropolis point 

is located at the intersection of the reversibility line with the first bisectrix (cr = (5). 
For 7 = 4/5, as in figure 11 we thus get 6 = —2/7. At this point, unfortunately, the 
quantities whose decay rates can be calculated by the above reduction technique are 
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the staggered magnetization and the staggered three-spin product [40], which are not 
relevant to the dynamics of local observables in a ferromagnetic model. 

To close, let us notice that the quantum Hamiltonian (5.2) does not have any 
special feature at the KDH and Metropolis points. In particular it seems to remain 
non-integrable. Further considerations on the subtle relationship between partial 
solvability and integrability in quantum chains and related matters can be found 
in [41] . The explicit examples analyzed there correspond to cases of free fermions. 
They are thus analogous to our solvable models. 


6. Discussion 

The present work contains a detailed investigation of the generic dynamics of the 
ferromagnetic Ising chain introduced in [1 [3 Ej. This generic one-dimensional 
kinetic Ising model extends the Glauber model [Tj to a two-parameter space 
corresponding to non-linearity and irreversibility, respectively measured by the 
deformation parameters 6 and g. While the introduction of the non-linearity 
parameter S is already present in [I] , the introduction of the irreversibility parameter g 
is a novelty. This generic dynamics is also the most general single-spin-flip 
dynamics which fulfills, besides global balance with respect to the ferromagnetic 
Hamiltonian spin reversal symmetry and a spatially homogeneous dynamical 

influence of nearest neighbours only. 

The present study has the virtue of organizing many partial results which were 
scattered in the literature concerning the role of the non-linearity parameter S for the 
reversible chain, and of unveiling novel features when the irreversibility parameter g 
is simultaneously present. In view of the extensive body of knowledge on kinetic Ising 
models, especially in one dimension, since Glauber’s seminal work in 1963 (see e.g. the 
reviews in [42] 1. it may appear surprising that only recently has this generic kinetic 
Ising model been considered. 

At infinite temperature, where Glauber dynamics accounts for independent spins, 
the presence of deformation parameters already has drastic consequences. The key 
observable is the overlap C{0, t) between the initial spin configuration and the current 
configuration at time t. Along the reversibility line, the non-linearity slows down the 
dynamics. The dependence of the relaxation rate of (7(0, t) on 5 is given exactly. An 
extreme non-linearity parameter (i5 = —1) yields a microcanonical dynamics, where the 
dynamics of domain walls is described by a SEP. The stretched exponential relaxation 
of (7(0, t) is also fully characterised. All along the reversibility line, an infinitesimal 
amount g of irreversibility induces an oscillatory relaxation of (7(0, t), with the period 
of oscillations diverging as 1/|??|. With irreversibility, the SEP domain wall dynamics 
becomes an ASEP dynamics, featuring a stretched exponential relaxation modulated 
by oscillations. Investigating the spectrum of the Markov matrices of finite chains 
directly yields a clear picture of the main features of the dynamics, such as reversibility 
or integrability (see figures [3 to [l^) ■ 

At finite temperature, the main novel feature is the occurrence of two successive 
thresholds in the irreversibility parameter g. This phenomenon was first uncovered by 
analytical means on the solvability line 5 = 0 [3] . Beyond a first threshold {g = ±r?c); 
a random and a thermalized initial states yield different spin relaxation times. This 
threshold also marks the onset of a strong violation of the equilibrium fluctuation- 
dissipation theorem. Beyond a second threshold {g = a-U two-time correlation 

functions exhibit an oscillatory relaxation. Figure [^summarizes the phase diagram in 
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the 7-77 plane for these correlation functions. An analogous phase diagram is expected 
to prevail for generic parameter values. In the present work we restricted the study 
to the second threshold, whose presence was demonstrated by simulations. 

The qualitative picture sketched above has been complemented by a good deal 
of quantitative predictions, coming from a variety of approaches including numerical 
simulations, time series expansions, and the spectra of Markov matrices and quantum 
spin Hamiltonians. 
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